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When looking for analytical approaches to treat frustrated quantum magnets, 

Qjl ' it is often very useful to start from a limit where the ground state is highly 

degenerate. This chapter discusses several ways of deriving effective Hamil- 

ti ' tonians around such limits, starting from standard degenerate perturbation 

theory and proceeding to modern approaches more appropriate for the deriva- 

C^ ' tion of high-order effective Hamiltonians, such as the perturbative continuous 

unitary transformations or contractor rcnormalization. In the course of this ex- 
position, a number of examples taken from the recent literature are discussed, 
including frustrated ladders and other dimer-based Heisenberg models in a 

Q ' field, as well as the mapping between frustrated Ising models in a transverse 

O . field and quantum dimer models. 
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As emphasized several times throughout this book, frustrated magnets often 
have a highly degenerate ground state in the classical limit. This is sometimes 
even taken as a definition of highly frustrated magnets. This degeneracy makes 
the semi-classical expansion in 1/S effectively uncontrolled (if it does not al- 

^^ ' ready fail simply because of divergent quantum fluctuations) , because usually 

one cannot perform and thus compare the expansions around all classical 
ground states. An infinite degeneracy is also often present in other limiting 
cases such as decoupled local units (such as triangles in the S ~ 1/2 trimerized 

^^ ' kagome lattice) or the Ising limit (for systems such as the antiferromagnctic 

Heisenberg model on the triangular or the kagome lattice). In such limits, 
which preserve the quantum nature of the problem, this degeneracy is not 
the end of the story, but rather the starting point of a systematic expansion, 
namely degenerate perturbation theory, which leads to an effective Hamilto- 
nian. In the context of strongly correlated systems, this type of method usually 
goes by the name 'strong-coupling expansion,' because the starting point of 
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the pcrturbative expansion is a Hamiltonian where only the interaction terms 
are kept, the kinetic terms being treated as the perturbation. 

There are several ways to perform this expansion, or more generally to 
derive an effective Hamiltonian. There are in fact two types of effective Hamil- 
tonian: those which act only in the degenerate subspace of a non-perturbed 
Hamiltonian, and those which act in the full Hilbert space but, through a 
canonical transformation, are rewritten as a series in terms of the ratio of two 
parameters. While the first type can always be written down explicitly, the 
second type can be derived in a simple way only provided that one can find 
a suitable generator. In the following, we will discuss both types of effective 
Hamiltonian, starting with the expansion in the degenerate subspace because 
this is more standard. 

The derivation of an effective Hamiltonian is extremely useful in isolating 
the relevant degrees of freedom. However, the problem is usually not solved 
once the effective Hamiltonian has been derived. Indeed, the new Hamiltonian 
often poses a problem as difficult as the original one. The primary advantage of 
the effective Hamiltonian is that, because the relevant degrees of freedom have 
been selected, simple approximations to the problem defined by the effective 
Hamiltonian often give deep insight into the physics of the problem. We will 
illustrate this point with several examples. 

This chapter is organized as follows. In Sec. II, we review briefly the degen- 
erate perturbation theory approach to effective Hamiltonians, with a concise 
but self-contained discussion of the second-order result, and a description of 
the form the expansion takes at higher orders. In Sec. HI, we discuss three 
examples taken from the field of frustrated magnetism where this approach 
has proven very useful: coupled dimers in a magnetic field, the Ising model 
in a transverse field, and the trimerized, spin-1/2 kagome antiferomagnet. In 
Sec. IV, we review more sophisticated approaches based on the same founda- 
tion: canonical transformations, continuous unitary transformations (CUTs), 
and the contractor renormalization group approach (CORE). Note that linked 
cluster expansions are discussed in another chapter of this book (that by 
Lauchli). We conclude in Sec. V with a discussion that includes a comparison 
of the various approaches. 



1.2 Strong-coupling expansion 

Let us consider a system described by a Hamiltonian 

acting in a Hilbert space H such that the ground state of Hq is degenerate. 
We denote by Hq the Hilbert space of the ground-state manifold. The goal is 
to find an effective Hamiltonian Hcs acting in Hq such that 

H,s\^)=E\^)^H\^)=E\^), \cf>)eHo, \4')eH. 
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1.2.1 Second-order perturbation theory 

Up to second order, the relation to be derived below is a standard result of 
quantum mechanics [1] . Denoting by Eq the ground-state energy of Hq and 
by Em the other eigenenergies, Ho\m) = Em\m) for Em ^ Eq. Thus for two 
vectors {(p), \(j)') G Ho, and up to second order in V, 

^ Eq - Em 

This result can be reformulated as an operator identity [2] . If one denotes by 
P the projector on 'Ho, and defines Q = 1 — P, then to second order in V 

PH,sP = PHoP + PVP + PVQ- 1-—-QVP 

Proof: suppose H\ip) = E\ip). Because P + Q ^ 1. this can be written as 

{P + Q)H{P + QM=E\^). 
Projecting onto Hq and V. — Ho gives 

PHP\i,)+PHQ\i))=EP\iP), (1) 

QHP\ij) + QHQli') = EQ\i,), (2) 

whence 

(2) => Q\^j) = {E- QHQ)-'QHP\^), 

(1) ^ PHP\i') + PHQ^-^j^QHP\i,) = EP\^). 

Expansion of {E - QHQ)-^ using {A - B)-^ = A-^ Y.n=o (-8^^^)": with 
A^ Eo- QHoQ and B = QVQ - E + Eq, leads to 

Pi/cffP = PHoP + PVP + PVQ{Eo - QHoQ)-^ 

OO 

E {{QVQ -E + Eo){Eo - QHoQY^YQVP. (1.1) 



OO 

X 

n=0 



The expansion is then truncated at n = 0. QED. 

1.2.2 High-order perturbation theory 

This form of the expansion is not well suited to the derivation of higher- 
order expansions, because beyond second order it contains explicitly the exact 
eigenenergy E. An expansion only in terms of the unperturbed eigenenergies 
can nevertheless be derived. The first systematic method for this dates to the 
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work of Kato [3]. Here wc follow the formulation of Takahashi [4], in which 
the expansion takes the form 

r = pp{ppp)-^/^, 

^PPPyi/2 = P + f2 i?!y^[F(P - P)Pr, 



n=l 



(2n)!! 



P = P-J2 J2 S'''VS''"'V...VS'''^+\ 

n—1 ki -{-... -\-kn-^-i— 71, ki>0 

S" = -P, s' ^ ^ 



^Eq — QHoQ J 
The true eigenstates ^ are related to the eigenstates (f) of -ffeff by 

\4')=r\i^), 

and likewise the observables transform according to 

o -^ r^'or. 

Thus the rtth order term of Hcs has the form 

^^ff'= Y. f{ki,k2,...,kn-l)VS'''VS''W...S'''^-'V, 

fei+...+fe„_i=n-l,fci>0 

where the coefficients /(fci, . . . , fc„_i) are deduced by appropriate bookkeeping 
from the previous expansions. The number of terms in such a strong-coupling 
expansion grows exponentially with n. In practice, beyond the fourth order it 
can generally be carried out only with the help of a computer. An alternative 
formulation based on continuous unitary transformations, which is simpler 
when applicable, will be discussed in the next section. In the remainder of this 
section, we discuss a number of selected examples where low-order degenerate 
perturbation theory provides considerable additional insight into the problem. 

1.2.3 Examples 

In the field of quantum magnetism, the best known example is the derivation 
of the Heisenberg model starting from the half-filled Hubbard model. The 
Hubbard model is defined by 

H = V + Ho^^t Y^ {c\^Cj„ + h.c.) + UY "*t"4- 
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Fig. 1.1. Left: schematic representation of the spin ladder of Eq. 1.2. The couplings 
entering Eq. 1.2 are denoted by the different line types: thick solid (J±), thin solid 
(Jii), dashed (J22), dot-dashed (J12) and dotted (J21). Right; magnetization curve 
of a frustrated spin ladder with a 1/2-magnetization plateau. 



At half-filling, the ground state of the interaction term Hq is 2^-fold de- 
generate, where TV is the number of sites. Treating the kinetic term y as a 
perturbation leads at second order (up to a constant) to the effective Hamil- 
tonian 

with J = At'^/U [5]. As noted in the introduction, this is a case where the 
effective Hamiltonian is itself very difhcult to solve, which is indeed true for the 
Heisenberg model on most lattices. As we shall see in the following, it is often 
useful to go one step further, and to derive a further effective Hamiltonian 
starting from one physically relevant limit. 

Example 1: frustrated spin-1/2 ladder in a magnetic field 

We consider the Heisenberg model for a frustrated spin-1/2 ladder in a mag- 
netic field defined by the Hamiltonian 

H =2_^J± Sin ■ S2n - B 2_^(5i„ + S2n) + /_, /_, ^ij ^ in ' Sjn+l- (1.2) 
n n n iJ — 1,2 

In the spin operators Sin, the index i refers to the leg and the index n to the 
rung (Fig. 1.1). The goal is to derive an effective Hamiltonian that describes 
the magnetization process in the limit J± ^ Jij [6,7]. 

The starting point of the perturbative expansion is the Hamiltonian of 
isolated dimers at the critical field Be = J± where, for one dimer, the triplet 
polarized in the field direction crosses the singlet, 



Ho — 2_^ J±Sln ■ S2n — Be / X^ln + ^^ 



(1.3) 



6 Frederic Mila^ and Kai Phillip Schmidt^ 

The ground state of this Hamiltonian is 2^-fold degenerate, where N is the 
number of ladder rungs, and the full Hamiltonian ean be treated within de- 
generate perturbation theory, the perturbation being given by 

^ - E E J^^ ^» ■ ^^™ - (^ - B,)Y.(Sin + SD- (1-4) 

(run) ' J = l,2 n 

The sum over (nm) refers to the nearest-neighbor rungs. A perturbative ex- 
pansion can be performed under the condition that the matrix elements of V 
are small compared to the excited states of Hq, i.e. as long as J^-, {B — B^) ^ 
J±. The condition on B — Be might suggest that such a calculation cannot 
give access to the full magnetization curve, but in fact the magnetization is 
rigorously equal to zero or to the full saturation value outside a window whose 
width is of order Jij . Thus this type of perturbation theory can indeed give 
the full magnetization curve. 

To write down the effective Hamiltonian, one needs a description of the 
Hilbert space. Because there are two states per rung, one simple choice is to 
introduce the Pauli matrices cr„ at each rung such that the singlet corresponds 
to \a^ = —1/2) and the triplet polarized along the field to \a^ = 1/2). In this 
basis, and up to first order in V, the effective Hamiltonian is given by 

i?cff - r' E(«+i + «+^^ + ^' E «+i - ^^ff E«) + ^' (1-5) 

n n n 

with 

J"^ = ^(Jll+J22-Jl2-J2l), (1.6) 

^■" = ^(Jll+J22 + Jl2 + J2l), (1.7) 

B,s = B-B,- i(,/n + J22 + ./12 + J21), (1.8) 

C = liJii + J22 + J12 + J21) -{B- B,). (1.9) 

o 

The effective Hamiltonian is identical to that of an XXZ chain in a field, which 
is a major step forward with respect to the original problem: this model has 
been investigated at length using the Bethe ansatz [8] and by field-theory 
methods [9, 10], and much is known about its low-energy properties. In par- 
ticular, there is a quantum phase transition in zero field at J^ = J^^ between 
a gaplcss phase at small J^ and a gapped phase at large J^ . 

To understand the physics of this phase transition for the original problem, 
it is expedient to perform a Jordan- Wigner transformation [11] of the effective 
model. In this language, the Hilbert space is that of spinless fermions on a 
chain, and the elementary operators are fermion creation (cj), annihilation 
(ci) and density (n^ = c\ci) operators at site i. An empty site corresponds to 
a rung singlet, an occupied one to a rung triplet, and the effective Hamiltonian 
becomes 
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^eff = -t^(c|cj+i + clj^^Ci) + V^TliTliJ^l ~ fl^Tli, (1.10) 

i in 

with t — J^y 12. V = J^ , and ji = B^a + J^ ■ In this model, the gapped phase 
of the ladder is a half-filled insulating phase of the fermionic chain, while the 
gapless one is a metallic phase (a Luttingcr liquid in this one-dimensional 
(ID) system). Thus the density as a function of the chemical potential has 
a plateau, the width of which is equal to the gap of the insulating phase. In 
the original model, this implies that the magnetization can have a plateau 
for certain parameters, a result which has been confirmed by density-matrix 
renormalization-group (DMRG) calculations [12]. The physics of the plateau 
state is discussed elsewhere in this volume (the chapter by Takigawa and 
Mila). For the purposes of the present chapter, we note only how powerful the 
effectivc-Hamiltonian method can be: a very simple, first-order calculation 
can basically solve the problem by mapping it onto another non-trivial but 
well-understood one. 

Example 2: expansion around the Ising limit 

In the previous example, as in the case of the Hubbard model, the unperturbed 
Hamiltonian is a sum of local terms, and the macroscopic ground-state de- 
generacy is given simply from the ground-state degeneracy of each term. This 
is not the only case where the effective-Hamiltonian approach is useful. An- 
other important example is the ground state of the antiferromagnetic Ising 
model, which is degenerate on non-bipartite lattices such as the triangular 
and kagome geometries. Again this can be the starting point of a degenerate 
perturbation theory towards the Heisenbcrg model on the same lattice if the 
transverse exchange is treated as a perturbation. In the same spirit, degenerate 
perturbation theory can be used to treat the effect of a transverse magnetic 
field applied to a frustrated Ising model. In this section we concentrate on the 
Ising case, which is of direct relevance to Quantum Dimcr Models (QDMs, 
discussed in the chapter by Moessner and Raman). 

Consider the fully frustrated Ising model in a transverse field on the hon- 
eycomb lattice, defined by the Hamiltonian 

H = Ho + V = -JY, M.ja^a] - T ^ < . (1.11) 

ihJ) ' 

In this expression, the parameters My = ±1 are chosen in such a way that 
their product around each hexagonal plaquette is equal to —1. All models 
in this class are then equivalent up to a gauge transformation, af — >■ eicrf, 
Mij -^ eiejMij, with a = ±1. One possible choice is shown in Fig. 1.2. 
For this particular choice, the two ferromagnetic states, with all spins either 
up or down, are ground states of Hq. Indeed, they satisfy five bonds out 
of every six, and it is clearly impossible to do better for a case with one 
antiferromagnetic coupling out of six. From a ferromagnetic ground state, one 
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Fig. 1.2. Left: representation of the model of Eq. 1.11. The thick vertical lines on 
the honeycomb lattice correspond to Mij = —1, all others to Mij — 1. The thick 
horizontal lines on the triangular lattice correspond to dimers in the ferromagnetic 
ground state. Right: dimer covering obtained from the ferromagnetic ground state 
after flipping the two spins shown as black dots. 



may then construct more ground states, and it is easy to verify that flipping 
the spins at the ends of a satisfied bond that connects two unsatisfied bonds 
leads to another ground state (right panel, Fig. 1.2). 

The structure of the Hilbert space is best understood by considering the 
dual lattice, which is the triangular lattice of sites at the centers of the hexag- 
onal plaquettes, and to draw a line between two neighboring sites of the dual 
lattice if the bond (of the original lattice) which it crosses is unsatisfied. If 
one imposes the constraint that each site be connected to exactly one of its 
neighbors, this set of lines defines a dimer covering of the triangular lattice, 
and there is a one-to-one correspondence between these dimer coverings and 
the ground states of Hq which have a given spin orientation for each given 
site. Because there are two possible spin orientations, the degenerate mani- 
fold of Hq is twice as large as the Hilbert space of dimer coverings, the two 
configurations leading to the same dimer covering being related by a global 
flip of all spins. 

For simplicity, we focus on the problem defined by the Ising Hamiltonian 
acting in the Hilbert space where two configurations related by a global spin- 
flip are identified. Then the degenerate subspacc can be described by the 
set of dimer coverings, and the effective Hamiltonian obtained by degenerate 
perturbation theory takes the form of a QDM. To deduce the form of this 
QDM, we first note that the effective Hamiltonian vanishes to first order in 
the transverse field. Indeed, fiipping a single spin changes the nature of the 
three bonds connected to it, which increases the number of unsatisfied bonds 
by one or by three, depending on whether or not one of the bonds connected 
to this site was unsatisfied. To second order, the situation is still more complex 
unless the two sites are nearest neighbors, in which case the process does not 
result in a state outside the ground-state manifold provided the two sites are 
on a satisfied bond connected to two unsatisfied ones. The resulting process 
leads to a dimer flip around a square plaquette of the dual lattice, and the 
effective Hamiltonian is a special case of the Rokhsar-Kivclson model, defined 
by the Hamiltonian 
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«E(|//>(//| 



H.c. 



(1.12) 



with t = r'^/J and w == 0. 

In this example, it is the effective model which is of direct interest in 
the context of frustrated magnets (as discussed in the chapter by Moessner 
and Raman) . The connection with the Ising model turns out to be helpful in 
identifying a possible phase transition, as first noted in Ref. [13], and led more 
recently to an analytical description of the fractional excitations which exist 
in the Resonating- Valence-Bond phase of the QDM on the square lattice [14]. 

As mentioned at the beginning of this section, an expansion around the 
Ising limit can also be performed for the XXZ version of the Heisenberg model, 
using the ratio J'^'^ j J^ as a small parameter. This approach has been used 
for the 1/3-plateau of the spin-1/2 kagomc antiferromagnet [15], for the half- 
magnetization plateau phase of the spin-1/2 pyrochlore Heisenberg antiferro- 
magnet [16], and for the 1/3-plateau phase of molecular analogs of the spin 
1/2 kagome antiferromagnet [17]. 



Example 3: weakly coupled triangles 

Frustration naturally appears in antiferromagnets in which the exchange paths 
create loops of odd length, the simplest of these loops being the triangle. 
There are several types of lattice geometry which can be considered as coupled 
triangles, and for which a perturbation theory starting from non-interacting 
triangles has proven to be useful. The effective model takes a special form due 
to the pecularity of the ground-state manifold of a triangle, which is four-fold 
degenerate. This is easily seen by rewriting the Hamiltonian in terms of the 
total spin of the triangle, ^tot = Si + -S'2 + ^3, to obtain 



H = J{Si • ^2 + ^2 • ^3 



s.i = i 



(Stot)' 



(1.13) 



a result emerging because it is possible to construct two doublets, i.e. four 
states, using three spin-1/2 entities. A convenient basis is provided by the 
simultaneous eigenstates of the scalar chirality, Si.{S2 x S3), and of the 
projection S^^^. of the total spin. 



\R, (t) = (I — aaa) + a;|a" — aa 
,2 



UJ \<T(T 



a»/x/3, 
L, ct) = (I — (Tcrcr) + w^lcr — aa) + uj\cra — cr))/v3. 



where ut = cxp(2i7r/3), cr = ±1/2 refers to S^^^., and L and R represent left- 
and right-handed chirality. 

In a system of weakly coupled triangles, treating the inter-triangle coupling 
as a perturbation leads to an effective Hamiltonian which acts in a Hilbert 
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space of dimension 4^', where Nt is the number of triangles. For a given 
triangle i, it is convenient to introduce a spin- 1/2 operator cr^, acting on the 
total spin, and a Pauli-matrix vector r^ acting in chirality space. To first order 
in the inter-triangle couplings, provided these couplings are SU(2)-invariant, 
the effective Hamiltonian then takes the general form 



E''^^-'^^^5' (1-14) 



HeS — Jeff 



where the sum is over the lattice of sites representing the arrangement of the 
triangles, Jeff is linear in the inter-triangle couplings, and the operator H^, 
depends on which sites of the triangles i and j are coupled by the inter-triangle 
interaction. 

In the case of a three-leg 'spin tube,' the effective model takes the explicit 
form [18] 

Hoff = y5]'<T, . a,[l + 2(t[t^, + r^r^J], (1.15) 

where J' is the rung coupling of the real model. For the effective model, 
field-theory arguments based on bosonization show that the spectrum must 
be gapped in all sectors, a non-trivial prediction to be contrasted with the 
gapless spectrum of the three- leg ladder [18] . 

Another example where this type of effective Hamiltonian has provided 
additional valuable insight is on the kagome lattice. Considering the trimer- 
ized version of this model, where the exchange constant within up-pointing 
triangles is taken to be J and in down-pointing triangles J', the effective 
Hamiltonian in the limit J' <^ J can be expressed as [19,20] 

J' ' 

iv) 

where the vectors e^ are taken from among ei = (1, 0), 62 = (^1/2, — \/3/2), 
and 63 = (— 1/2, •\/3/2), the choice depending on the labelling of the spins 
inside each triangle [21]. For this effective Hamiltonian, a mean- field decou- 
pling of spin and chirality leads to a highly degenerate ground state. Each 
ground state can be associated with a dimer covering of the triangular lat- 
tice of up- (or down-)pointing triangles, and the number of states accordingly 
grows as 1.5351^' in terms of the number of triangles, or as 1.1536^ in terms 
of the number of spins [20] . This compares well with the number of low- lying 
singlets observed in numerical simulations of the spin-1/2 kagome antiferro- 
magnet [22]. 

This effective Hamiltonian is also a very useful starting point to discuss 
the physics in a magnetic field. Antiferromagnets composed of weakly coupled 
triangles exhibit a plateau at magnetization 1/3, in which all triangles have 
a total spin equal to 1/2 (in a S" = 1/2 system) and oriented in the field 
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direction. Inside the plateau, cri ■ ctj ~ 1/4 is a constant, and the effective 
Hamiltonian is a pure chirality model. This plateau has been investigated for 
the kagome lattice in Ref. [23]. 

Finally, similar ideas based on weakly coupled tetrahedra have provided 
equally helpful insight into the properties of the spin-1/2 Heisenberg model 
on the pyrochlore lattice and in related systems [24,25]. 



1.3 Alternative approaches yielding effective 
Hamiltonians 

There are several ways of deriving effective Hamiltonians based on techniques 
other than direct strong-coupling expansions. The aim of this section is to 
provide a review of the physics underlying these approaches, with appropriate 
references for further reading concerning their more detailed implementations. 

1.3.1 Canonical transformation 

The canonical transformation of a Hamiltonian is defined by iJ — > i7 = 
e^He~^ , where 77 is antihermitian, so that e^' is unitary. If [^) is an eigenstate 
of H, then \ip) = e''['0) is an eigenstate of H with the same eigenvalue. If 
the operators are transformed simultaneously according io A ^f A = e^ Ae~'^ , 
then (^[i[^) = {ip\A\i;). 

The foundation for using a canonical transformation to derive an effective 
Hamiltonian is the identity 

e"ife-'' = F+[r/,ff] + l[r?,[r?,/f]] + ... 



00 ^ 



n 

n=l 



which is simply the Taylor expansion of H{X) = e^''He^^^' for A = 1. 

Considering the case where H = Hq + XV , if one can find an operator 77 
such that [?7, Hq] = —V, then using A77 as a generator leads to 

This operator is a series in powers of A, and hence of the perturbation V . 
While its structure is reminiscent of the results of high-order perturbation 
theory, there is an important difference: H acts in the full Hilbert space of 
H, whereas H^s acts in the ground-state manifold of Hq. Depending on the 
problem, this may or may not be an advantage. If one is interested in high- 
energy states which may be detected in a particular experiment, the canonical- 
transformation approach has distinct advantages, because it gives all of the 
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eigenstatcs up to a certain order, and not only those which have evolved from 
the ground-state manifold of H^ under the perturbation. By contrast, if the 
goal is to reduce the Hilbert space to study the low-energy sector using, for ex- 
ample, exact-diagonalization calculations on finite clusters, then a degenerate 
perturbation theory is sufficient. 

1.3.2 CUT (Continuous Unitary Transformation) 

The canonical transformation introduced in the previous section is by no 
means the only possibility for obtaining an effective model by a unitary trans- 
formation. In fact there are many ways to do this, even in low orders of a 
perturbative approach, and it is therefore an obvious question to ask whether 
an optimal choice of transformation exists. 

This question led both Wegner [30] and Glazek and Wilson [31,32] to in- 
troduce independently of each other the concept of continuous unitary trans- 
formations (CUTs) [33]. In contrast to the one-step transformation discussed 
in the last section, here the unitary transformation is constructed as an infi- 
nite product of infinitesimal transformations. Although measurable (on-shell) 
quantities, such as energies, have to be the same independent of which kind 
of transformation has been chosen, off-shell quantities such as effective inter- 
actions can differ strongly. This has been demonstrated in a quite impressive 
manner for the case of the Frohlich Hamiltonian, which describes conventional 
superconductivity mediated by the electron-phonon interaction [34] . Here the 
effective electron-electron interaction at second order in the electron-phonon 
coupling shows divergences for the case of a one-step transformation (previ- 
ous section), whereas in the continuous treatment the attractive interaction 
is smooth. 

Another respect in which the approaches differ is that the one-step canoni- 
cal transformation can be applied practically only at low orders in the pertur- 
bation. However, there are physical situations where one is interested in the 
quantitative determination of an effective Hamiltonian for a given parameter 
set in the original model. One example is the case of strongly frustrated net- 
works of coupled dimers, as in the Shastry-Sutherland model, where processes 
relevant to the magnetization of the system appear only at high orders. 

In the method of CUTs, a continuous parameter / is introduced such that 
I = refers to the initial system H and / = cx3 corresponds to the final 
effective system, which should correspond to a simplified physical picture. The 
transformation can be constructed such that processes at higher energies are 
treated before those at lower energies. This renormalizing property is similar 
to Wilson's renormalization-group approach [35]. 

Let U be the unitary transformation which diagonalizes the Hamiltonian 
H and let H{1) ~ U^{l)HU{l). This unitary transformation is then equivalent 
to performing an infinite sequence of unitary transformations, e"'''-'-*'*', with 
the antihermitian generators 
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r^il) ^ -UHl)diU{l). (1.17) 

Taking the derivative with respect to / results in the "flow equation" 

diH{l)^[iil),H{l)], (1.18) 

which defines the change of the Hamihonian during the flow. Note that 
Eq. 1.18 represents an infinite hierarchy of coupled differential equations, be- 
cause for the general case an increasing number of terms is generated on the 
right-hand side at each order. In practice, one has therefore to perform a 
truncation (below). 

The properties of the effective Hamiltonian depend strongly on the choice 
of the generator. There are in essence two different modern variants. The first 
one uses the generator introduced originally by Wegner [30], which aims to 
eliminate interaction matrix elements with the goal of obtaining an energy- 
diagonal effective Hamiltonian. This approach has been applied successfully 
to a large class of problems, with special attention being given to determining 
the ground state of interacting quantum many-body problems [30,34,36]. 
The second variant is the quasiparticle-conserving CUT which, as its name 
suggests, maps the Hamiltonian Hq to an effective Hamiltonian conserving the 
number of quasiparticles [37-40]. This approach can be used either to study 
the excitations of an already known quantum ground state [33,40-44], one 
application for which has been to bound states, or, in analogy to the previous 
case, to derive effective low-energy models [45-49]. 

Returning to the point mentioned above, in all practical calculations it is 
necessary to truncate the flow equation (1.18). For this there are two options: 
i) cutting the hierarchy at one level and solving the equations numerically 
(self-similar CUTs) or ii) using a series-expansion ansatz for rj and H, and 
solving the flow equations perturbatively to high order (perturbative CUTs). 

Here we focus only on presenting one illustrative example, for which we 
choose the perturbative version of quasiparticle-conserving CUTs [37,40,50, 
51]. If the problem at hand meets the two conditions: 

1. the unperturbed part has an equidistant spectrum bounded from below, 

2. there is an integer number iV > such that the perturbing part can be 
separated as X]n=-Ar^n7 where r„ increments (or if n < decrements) 
the number of particles by n, 

then the CUT in its quasiparticle-conserving form can be solved to high order 
in the perturbation and the effective Hamiltonian is given by the general 
form [40] 

oo 

H,s{x) = Q + "^x'' Yl <^in^)T{m), (1.19) 

fc=l ImNfc 



where Q is the unperturbed part of the Hamiltonian, x an expansion param- 
eter, m = (toi, ?7i2, . . . , ?7Zfc), and M(m) = X]i=i "^i ~ ^ reflects the con- 
servation of the number of particles. The action of Heft can be viewed as a 
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weighted sum of particle-number-conscrving virtual excitation processes, each 
of which is encoded in a monomial T{rn) = Tm^Tm^ . . -Tm^. The coefficients 
C{rn) are rational numbers which can be calculated (to high order in the 
perturbation) exactly as the ratio of two integers. It should be emphasized 
that the effective Hamiltonian HcS, which has known coefficients C{Tn), can 
be used straightforwardly in all perturbative problems that meet the above 
conditions. 

We illustrate the method first for an unfrustrated spin ladder in a magnetic 
field, including from Eq. (1.2) only the magnetic exchange J± on the rungs, 
the magnetic field B, and the unfrustrated exchange along the legs of the 
ladder, setting Jn = J22 = Jii. Rewriting the Hamiltonian in terms of rung 
triplet operators i^, with a ~ {x, y, z} [52], we obtain 

H = J±Q + ^ [To + T+2 + T_2] + Hb, 
where 

i,a 

i,a 
-^—2 / ^ ^a.i^a.i-\-l — +2" 

The operator Q counts the total number of triplet excitations, while the oper- 
ators r„ change the triplet number by n, and Hb denotes the magnetic-field 
term. 

In the following we consider the limit of weakly coupled rung dimers, i.e. 
we set Jx = 1 and consider Ju/J± = x as the small expansion parameter. 
The effective Hamiltonian i/cff obtained by a quasiparticlc-conserving CUT 
has the property [_ffoff,Q] = 0, meaning that the total number of triplons 
(dressed triplets which are the elementary excitations of the spin ladder) is a 
conserved quantity. The effective Hamiltonian at second order is 

2 

Uf^ = Q + xTo + ^ [r+2T_2 - ^-2^2] + Hb. 

The total spin S^^^ is a conserved quantity. The magnetic-field term is there- 
fore not changed under the unitary transformation, and the low-energy physics 
is influenced solely by the local singlet \s) and the triplet \t^) polarized par- 
allel to the magnetic field (as discussed in Sec. 2). Identifying \s) with an 
empty site and |t^) with the presence of a hard-core boson, one may deduce 
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the effective Hamiltonian in this basis by calculating matrix elements on a 
finite cluster [50,51]. The result is 

»je{i,2} i 

-t[ ^ f 6i_i&i+i + h.c.j Hi + Vj ^ n^Ui+j, (1-20) 

i ij6{l,2} 

with ti = xl2, t2 = x^/A, t[ = -a-V4, wi = x/2 - Sx^/g, V2 = 0, and 
11 = B - 1 + 3x74. 

We emphasize that a calculation such as this, based on coupled dimers in 
an external magnetic field, can be performed for any frustrated lattice and 
to high orders in the perturbation. By using appropriate extrapolations, a 
quantitative low-energy effective Hamiltonian may be derived, which is usually 
applicable over a large part of the parameter space. The high-order expansion 
(including extrapolation of the bare series) becomes problematic only when 
the correlation length of the system exceeds the spatial range covered by the 
maximum order treated, for example close to a quantum phase transition. The 
method will break down if the ground state for the parameters of interest is 
not unitarily connected to the ground state about which one is expanding. 

As a second example illustrating the importance of a quantitative effective 
model, meaning one obtained with high-order accuracy, we discuss the 2D 
spin-1/2 Heisenberg system known as the Shastry-Sutherland model [53] in a 
magnetic field, 

H = fY,S^■ S,+J Y. S^- S, -BY^St. (1.21) 

The bonds denoted <C i,j ^ represent an array of orthogonal dimers, 
while the bonds {i,j), which are intcr-dimcr couplings, form a square lattice 
(Fig. 1.3). This theoretical model is believed to be realized experimentally in 
the layered copper oxide material SrCu2 (603)2, where the coupling ratio is 
J' / J Ri 0.63. In the theoretical model, for J' /J smaller than a critical ratio 
of order 0.7, the ground state of the model is given exactly by the product of 
dimer singlets, and the magnetization process of the system can be described 
in terms of hard-core bosons which, as discussed for the spin ladder above, 
represent polarized \t^) triplons on the dimers, interacting and moving on an 
effective square lattice [57,58]. 

A consequence of the strong frustration is the appearance of several magne- 
tization plateaus which correspond to Mott-insulating phases of the hard-core 
bosons [54-56], where the translational symmetry of the system is broken 
and triplon excitations are frozen in the ground state as in a charge-ordered 
state (discussed in the chapter by Takigawa and Mila) . Theoretically, all ap- 
proaches to the basic theoretical model agree on the presence of magnetization 
plateaus at 1/3 and 1/2 of the saturation value [57-61], in agreement with ex- 
periments [55,62]. However, the structure below 1/3 magnetization is rather 
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controversial. On the experimental side, the original pulsed-field data show 
only two anomalies which were interpreted as plateaus at 1/8 and 1/4 [55], 
but the presence of additional phase transitions, and of a broken translational 
symmetry above the 1/8-plateau have been established by recent torque and 
NMR measurements performed up to 31 T [63,64]. The possibility of addi- 
tional plateaus has been pointed out by Sebastian et al. [62] , who interpreted 
their high-field torque measurements as evidence for plateaus at 1/g, with 
2 < g < 9, and at 2/9. On the theoretical side, the situation is also not 
settled. The finite clusters available for exact-diagonalization studies are not 
large enough to allow reliable predictions for high-commensurability plateaus, 
and the accuracy of the Chern-Simons mean-field approach initiated by Mis- 
guich et al. [60], and employed recently by Sebastian et al. [62] to explain 
their apparent additional plateaus, is not easy to assess. The essential diffi- 
culty lies in the fact that, because plateaus are a consequence of repulsive 
interactions between triplons, an accurate determination of the low-density, 
high-commensurability plateaus requires a precise knowledge of the long-range 
part of the interaction. 
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Fig. 1.3. Shastry-Sutherland lattice and definition of the two-body interactions. V„ 
is the coeflicient of the two-body interactions between the reference dimer (thick 
bond) and the corresponding dimer. Figure courtesy of Ref. [65]. 



Such a precise analysis was conducted recently using perturbative CUTs 
[65] . The processes relevant for the physics in a finite magnetic field are those 
with maximum total spin and total Sz • Thus the general form of the effective 
Hamiltonian obtained by the perturbative CUT takes the form 



HeS — 



j: E <^r 



A 



l}^ h 



(1.22) 



n=2,4,6... ri, 



where r, labels the sites of the square lattice formed by the J bonds, while 
the hard-core boson operator &J creates a polarized \t^) triplon at site r. 
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The coefficients Cri^r2,...,r„ a-rc then obtained as iiigh-order scries in J'/ J, 
computed up to order 15 for the two-body interactions Vn to be discussed 
below. 




Fig. 1.4. Coefficients of the extrapolated two-body interactions Vn as function of 
J' /J. Inset: different extrapolations (solid lines) shown together with the bare series 
(dashed lines) for V3 and V5 . Figure courtesy of Ref . [65] . 



It is found that the magnitudes of all the interaction terms decrease when 
the separation of the sites is increased. In addition, the physics at low density 
is dominated by the two-body density-density interactions, while the standard 
two-site hopping is, as expected, strongly suppressed due to the frustration 
[66]. The evolution with J' /J of the most relevant two-body interactions, 
defined in Fig. 1.3, is shown in Fig. 1.4. At small J' /J, interactions beyond V4 
are small and may be neglected, but for larger coupling ratios the higher-order 
terms V^, V5, and Vj (appearing at order 6) become important and contribute 
to the formation of low-density plateaus. 

In general, the effective Hamiltonian HcS is by no means simpler than 
the original one, but it becomes so in the limit of low density and moderate 
J'/ J. Indeed, in this limit the kinetic terms are very small, and they can 
be considered as a perturbation of the interaction part, which is diagonal in 
the local Fock basis, \nr-^ , Ur^ , ■ ■ ■) ■ It is thus appropriate to use a Hartree 
approximation in which the variational ground state is a product of local 
boson wave- functions, because this approximation becomes exact in the limit 
of vanishing kinetic energy. This approach has the further advantage that it 
can be used to compare and treat rather large unit cells, which is important in 
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the SrCu2 (603)2 problem, where "solid" phases with a complicated structure 
can be found at low magnetizations. 

The resulting phase diagram in the Hartree approximation is shown in 
Fig. 1.5. The phase diagram is dominated by a series of plateaus, which appear 
at 1/3 and 1/2 (not shown) even at very small J' / J , and at 2/9, 1/6, 1/9, and 
2/15 as J' / J is increased. A plateau structure of this kind is to be expected 
when the kinetic terms, as here, are very small, because if these terms were 
completely absent, the magnetization curve would consist simply of a sequence 
of plateaus. At J' / J = 0.5, the 1/6-plateau is by far the most prominent 
stucture below 1/3. 

All of the plateaus found at low densities are actually stabilized by two- 
body repulsive interactions Vn appearing at high orders in J' / J . It is therefore 
crucial to obtain the effective Hamiltonian very accurately, because the most 
significant features in this density regime are found to result from the compe- 
tition between small interactions. 

Finally, it is worth adding that the method of CUTs is also able to treat 
observables [36,50,65]. To this end, an observable O must be transformed by 
the same CUT as the Hamiltonian, 

diO{l) = W),0{l)], (1.23) 

yielding effective observables Ooff = 0{l = 00). Here we mention only one pos- 
sible application which can be very useful for the physics of frustrated quan- 
tum magnets: a typical situation in quantum magnetism is that, in addition 
to the dominant nearest-neighbor Heisenberg exchange interactions, there are 
small coupling terms ii/^add, such as Dzyloshinskii-Moriya interactions, which 
can have a profound influence on the physics of the system. Formally 

^tot = H + i/add , 

and it is both elegant and efficient to perform first a CUT on H, which contains 
the dominant couplings, and then to treat i^add as an observable 

Here H°^ conserves the number of quasiparticles, whereas the transformed 
observable U'^H^ddU does not. Thus in the second step, performed after the 
first CUT, either ordinary perturbation theory or a second CUT can be applied 
to treat the term [/^_ffaddC^, which contains the small couplings. In cases where 
-ffadd mixes low- and high-energy states, it is essential to retain access to the 
full Hilbert space of the problem. 

1.3.3 CORE (Contractor Renormalization) 

One further method which we highlight here for the derivation of low-energy 
effective models is the CORE technique, invented by Morningstar and Wein- 
stein [68] . The underlying idea is to derive effective Hamiltonians for a trun- 
cated local basis in such a manner that the low-energy spectrum of the model 
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Fig. 1.5. Magnetization plateaus as a function of ^ and J' I J . The boson density 
n is equal to the magnetization in units of the saturation value, and the chemical 
potential \x is equal to the magnetic field B. Solid lines denote results that are fully 
converged with respect to the terms retained in the Hamiltonian. Well-converged 
results are then connected by dashed lines. Figure courtesy of Ref. [65]. 



under study is reproduced exactly. As for CUTs, the CORE approach can be 
used either in an analytically oriented form [69-73] or as a numerical tech- 
nique [74-77]. For recent reviews on CORE we refer the reader to Rcfs. [78,79]. 
The essential steps of the CORE algorithm are the following. 

1) The system is divided into local subunits. One subunit is diagonalized, 
keeping M suitable low-energy states. 

2) The full Hamiltonian is diagonalized on a connected graph consisting of N 
subunits. The low-energy eigenenergies e„ and eigenstates |n) are calculated. 

3) A basis of dimension M^ is obtained by projecting the eigenstates onto 
the tensor product space of the retained states. 

4) The effective Hamiltonian is constructed accorchng to 



rroff 



M" 



i\i>n){'ipn 



5) The connected rangc-A^ interactions are determined by subtracting the 

contributions of all connected subclusters. 

Finally, the effective Hamiltionian is deduced by a cluster expansion as 



h: 



eff 
CORE 



= E^ + E^: 



(iJ) 



Note that ii^coRE reproduces exactly the low-energy physics if one considers 
all of the terms on the right-hand side. 

In practice, it is necessary to perform a truncation. The convergence of the 
algorithm therefore depends both on the range of the operators taken in the 
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cluster expansion and on the number and type of low-energy states retained 
for one subunit. Hence the successful application of the CORE technique does 
require some physical insight concerning the problem at hand. However, once 
the relevant degrees of freedom have been selected, CORE represents a non- 
perturbative method for deriving effective low-energy Hamiltonians. 

An important feature of the CORE algorithm is that it does not rely on the 
system being in a certain physical phase (to be contrasted with the example of 
quasiparticle-conserving CUTs discussed in the preceding section) and there- 
fore does not break down even if a quantum phase transition takes place in the 
parameter space of the original model. As an example, we present here some 
illustrative CORE results for the magnetization curves of the 2D spin-1/2 
Heisenberg model on the Shastry-Suthcrland lattice discussed in the previous 
section [67]. The CORE technique is used to derive an effective Hamiltonian 
which is then treated by exact diagonalization (ED). The effective model ob- 
tained by CORE is found to agree very well with that obtained by perturbative 
CUTs. It is therefore expected that differences between the two approaches 
arise primarily from the method used to treat the effective model, which in 
the examples shown is either a classical limit [65] or ED [67]. 

Results obtained from CORE-I-ED are shown in Fig. 1.6. As in the treat- 
ment by a perturbative CUT, a rich plateau structure is found below m = 1/3. 
This again highlights the utility of an approach to the physics by first deriv- 
ing a quantitative effective model which is then treated by simpler techniques. 
One obvious advantage of ED compared to the classical treatment is that it 
takes quantum fluctuations fully into account. It is therefore striking that the 
magnetization curves in both approaches are dominated by the presence of 
plateaus. By contrast, one drawback in using ED as a solver for the effective 
model (and one advantage of the classical solver) is the restriction on cluster 
sizes and shapes. Because the physics at low densities involves solid phases 
with rather large unit cells, it is a challenge within the ED approach to stabi- 
lize and to compare different solids, such as the 2/15 phase found in Ref. [65]. 
Despite the differing aspects, both positive and negative, of these approaches 
to this problem, it should be emphasized that the advanced techniques used 
to derive effective Hamiltonians have proven to be crucial in the discovery and 
resolution of the complicated magnetization processes of strongly frustrated 
quantum magnets. 



1.4 Conclusions 

In this chapter we have summarized a variety of tools which are used in the 
field of highly frustrated magnetism to derive effective low-energy Hamiltoni- 
ans. We have aimed to capture the essential technical aspects of these different 
approaches and to provide examples of them for a number of physical appli- 
cations. While it is not possible in a chapter of this type to cover all such 
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Fig. 1.6. (a) Magnetization curves obtained by ED and CORE calculations on an 
Af = 32 lattice and from CORE results with A/' = 64 for J' / J = 0.5. (b) Phase 
diagram for Ai" = 32 as a function of J' / J and magnetic field h/J: CORE results 
(lines) agree with ED (symbols) for locations of the m — 1/2, 1/4, and 1/8 plateaus 
which are allowed on this cluster. (c,d) Same as (a,b) for Af = 36 and A^ — 72. On 
these clusters, the m = 1/2, 1/3, 2/9, 1/6, and 1/9 plateaus are allowed. Figure 
courtesy of Ref. [67] . 



methods in full detail, where relvant we have referred the interested reader to 
the more specialized literature. 

We have shown that strong-coupling expansions and the derivation of ef- 
fective models are at the same time standard techniques used by theoretical 
physicists for a broad range of physical questions and also a very active area of 
current research on strongly correlated quantum systems. For the first type of 
application, it is often the goal to identify the relevant low-energy degrees of 
freedom and to use an effective model to leading order in perturbation theory 
to understand qualitatively the phase diagram of a given model. By contrast. 



22 Frederic Mila^ and Kai Phillip Schmidt^ 

the aim of the current developments is to obtain a quantitative derivation 
of effective models and a complete understanding of the breakdown of such 
derivations. 



References 

1. See for example L. Landau and E. Lifshitz, Quantum mechanics, (Pergamon, 
1991). 

2. See for example P. Fulde, Electron Correlations in Molecules and Solids, 
(Springer, 1991). 

3. T. Kato, Prog. Theor. Phys. 4, 514 (1949). 

4. M. Takahashi, J. Phys. C 10, 1289 (1977). 

5. P. W. Anderson, Phys. Rev. 115, 2 (1959). 

6. K. Totsuka, Phys. Rev. B 57, 3454 (1998). 

7. F. Mila, Eur. Phys. J. B 6, 201 (1998). 

8. H. A. Bethe, Z. Phys. 71, 205 (1931). 

9. A. Luther and I. Peschel, Phys. Rev. B 9, 2911 (1974). 

10. F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980). 

11. See for example D. C. Mattis, The theory of magnetism made simple, (World 
Scientific, 2006). 

12. J.-B. Fouet, F. Mila, D. Clarke, H. Youk, O. Tchernyshyov, P. Fendley, and R. 
M. Noack, Phys. Rev. B 73, 214405 (2006). 

13. R. Moessner, S. L. Sondhi, and P. Chandra, Phys. Rev. Lett. 84, 4457 (2000). 

14. G. Misguich and F. Mila, Phys. Rev. B 77, 134421 (2008). 

15. D. C. Cabra, M. D. Grynberg, P. C. Holdsworth, A. Honecker, P. Pujol, J. 
Richter, D. Schmalfuss, and J. Schulenburg, Phys. Rev. B 71, 144420 (2005). 

16. D. L. Bergman, R. Shindou, G. A. Fiete, and L. Balents, Phys. Rev. B 75, 
094403 (2007). 

17. L Rousochatzakis, A. M. Lauchh, and F. Mila, Phys. Rev. B 77, 094420 (2008). 

18. H. J. Schulz in Proceedings of the XXXIst Rencontres de Moriond, edited by 
T. Martin, G. Montambaux, and J. Trin Thanh Vin (Editions Frontieres, Gif- 
sur-Yvette, 1996). 

19. V. Subrahmanyam, Phys. Rev. B 52, 1133 (1995). 

20. F. Mila, Phys. Rev. Lett. 81, 2356 (1998). 

21. For a detailed discussion of this point, see M. Ferrero, F. Becca, and F. Mila, 
Phys. Rev. B 68, 214431 (2003). 

22. P. Lecheminant, B. Bernu, C. Lhuillier, L. Pierre, and P. Sindzingre, Phys. Rev. 
B 56, 2521 (1997). 

23. B. Damski, H.-U. Everts, A. Honecker, H. Fehrmann, L. Santos, and M. Lewen- 
stein, Phys. Rev. Lett. 95, 060403 (2005). 

24. H. Tsunetsugu, Phys. Rev. B 65, 024415 (2002). 

25. V. N. Kotov, M. E. Zhitomirsky, M. Elhajal, and F. Mila, Phys. Rev. B 70, 
214401 (2004); V. N. Kotov, M. Elhajal, M. E. Zhitomirsky, and F. Mila, Phys. 
Rev. B 72, 014421 (2005). 

26. H. J. Schmidt and Y. Kuramoto, Physica C 167, 263 (1990). 

27. M. Takahashi, J. Phys. C: Solid State Phys. 10, 1289 (1977). 

28. A. H. MacDonald, S. M. Girvin, and D. Yoshioka, Phys. Rev. B 41, 2565 (1990). 

29. E. Miiller-Hartmann and A. Reischl, Eur. Phys. J. B 28, 173 (2002). 



1 Strong-coupling expansion and effective liamiltonians 23 

30. F. J. Wegner, Ann. Pliysik 3, 77 (f994). 

3f. S. D. Glazek and K. G. Wilson, Pliys. Rev. D 48, 5863 (1993). 

32. S. D. Giazek and K. G. Wilson, Phys. Rev. D 49, 4214 (1994). 

33. C. Knetter, K. P. Schmidt, M. Griininger, and G. S. Uhrig, Phys. Rev. Lett. 
87, 167204 (2001). 

34. P. Lenz and F. Wegner, Nucl. Phys B 482, 693 (1996). 

35. K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

36. S. Kehrein, Flow-equation approach to many-particle systems, Springer Tracts 
in Modern Physics 217 (2006), and references therein. 

37. J. Stein, J. Stat. Phys. 88, 487 (1997). 

38. A. Mielke, Eur. Phys. J. B 5, 605 (1998). 

39. G. S. Uhrig and B. Normand, Phys. Rev. B 58, R14705 (1998). 

40. C. Knetter and G. S. Uhrig, Eur. Phys. J. B 13, 209 (2000). 

41. C. Knetter, A. Biihler, E. Miiller-Hartmann, and G. S. Uhrig, Phys. Rev. Lett. 
85, 3958 (2000). 

42. C. P. Heidbrink and G. S. Uhrig, Eur. Phys. J. B 30, 443 (2002). 

43. W. Brenig, Phys. Rev. B 67, 064402 (2003). 

44. K. P. Schmidt and G. S. Uhrig, Mod. Phys. Lett. B 19, 1179 (2005). 

45. A. Reischl, E. Miiller-Hartmann, and G. S. Uhrig, Phys. Rev. B 70, 245124 
(2004). 

46. S. Dusuel and G. S. Uhrig, J. Phys. A: Math, and Gen. 37, 9275 (2004). 

47. S. Dusuel and J. Vidal, Phys. Rev. Lett. 93, 237204 (2004). 

48. J. N. Kriel, A. Y. Morozov, and F. G. Scholtz, J. Phys. A: Math, and Gen. 38, 
205 (2005). 

49. K. P. Schmidt, S. Dusuel, and J. Vidal, Phys. Rev. Lett. 100, 057208 (2008). 

50. C. Knetter, K. P. Schmidt, and G. S. Uhrig, J. Phys.: Condens. Matter 36, 
7889 (2003). 

51. C. Knetter, K. P. Schmidt, and G. S. Uhrig, Eur. Phys. J. B 36, 525 (2004). 

52. S. Sachdev and R. N. Bhatt, Phys. Rev. B 41, 9323 (1990). 

53. B. S. Shastry and B. Sutherland, Physica B 108B, 1069 (1981). 

54. E. Miiller-Hartmann, R. R. P. Singh, C. Knetter, and G. S. Uhrig, Phys. Rev. 
Lett. 81, 1808 (2000). 

55. K. Onizuka et al, J. Phys. Soc. Jpn. 69, 1016 (2000). 

56. K. Kodama et al, Science 298, 395 (2002). 

57. T. Momoi and K. Totsuka, Phys. Rev. B 62, 15067 (2000). 

58. S. Miyahara and K. Ueda, J. Phys.: Condens. Matter 15, R327 (2003). 

59. S. Miyahara and K. Ueda, Phys. Rev. B 61, 3417 (2000). 

60. G. Misguich, T. Jolicoeur, and S. M. Girvin, Phys. Rev. Lett. 87, 097203 (2001). 

61. S. Miyahara, F. Becca and F. Mila, Phys. Rev. B 68, 024401 (2003). 

62. S. E. Sebastian, N. Harrison, P. Sengupta, C. D. Batista, S. Francoual, E. Palm, 
T. Murphy, H. A. Dabkowska, and B. D. Gaulin, Proc. Nat. Acad. Sci. 105, 
20157 (2008). 

63. F. Levy, I. Sheikin, C. Berthier, M. Horvatic, M. Takigawa, H. Kageyama, T. 
Waki, and Y. Ueda, Europhys. Lett. 81, 67004 (2008). 

64. M. Takigawa, S. Matsubara, M. Horvatic, C. Berthier, H. Kageyama, and Y. 
Ueda, Phys. Rev. Lett. 101, 037202 (2008). 

65. J. Dorier, K. P. Schmidt, and F. Mila, Phys. Rev. Lett. 101, 250402 (2008). 

66. S. Miyahara and K. Ueda, Phys. Rev. Lett. 82, 3701 (1999). 

67. A. Abendschein and S. Capponi, Phys. Rev. Lett. 101, 227201 (2008). 



24 Frederic Mila^ and Kai Phillip Schmidt^ 

68. C. J. Morningstar and M. Weinstein, Phys. Rev. Lett. 73, 1873 (1994); Phys. 
Rev. D 54, 4131 (1996). 

69. M. Weinstein, Phys. Rev. B 63, 174421 (2001). 

70. E. Altman and A. Auerbach, Phys. Rev. B 65, 104508 (2002). 

71. E. Berg, E. Altmann, and A. Auerbach, Phys. Rev. Lett. 90, 147204 (2003). 

72. R. Budnik and A. Auerbach, Phys. Rev. Lett. 93, 187205 (2004). 

73. M. S. Siu and M. Weinstein, Phys. Rev. B 75, 184403 (2007). 

74. J. Piekarewicz and J. R. Shepard, Phys. Rev. B 56, 5366 (1997); Phys. Rev. B 
57, 10260 (1998). 

75. S. Capponi and D. Poilblanc, Phys. Rev. B 66, 180503(R) (2002). 

76. S. Capponi, A. Lauchh, and M. Mambrini, Phys. Rev. B 70, 104424 (2004). 

77. A. Abendschein and S. Capponi, Phys. Rev. B 76, 064413 (2007). 

78. A. Auerbach, AIP Conf. Proc. 816 1 (2006). 

79. S. Capponi AIP Conf. Proc. 816 16 (2006). 



